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Abstract 

We report first principles density functional calculations of the spontaneous polarizations, piezo- 
electric stress constants and elastic constants for the II- VI wurtzite-structure oxides ZnO, MgO 
and CdO. Using our pseudopotential self-interaction corrected implementation of density functional 
theory, we obtain polarization values of -0.05, -0.17 and -0.10 C/m 2 , and piezoelectric constants, 
e 33 (e 3 i) of 1.34 (-0.57), 2.26 (-0.38) and 1.67 (-0.48) C/m 2 for structurally relaxed ZnO, MgO and 
CdO respectively. These properties are consistently larger in magnitude than the corresponding 
GaN, A1N and InN analogues. Therefore we predict that larger internal fields will be attainable in 
ZnO-based polarization field effect transistors than in equivalent GaN-based devices. 
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Over the past decade, nitride-based III-V semiconductors have emerged as leading con- 
tenders for many technological applications, particularly in optoelectronics (blue and green 
lasers and light emitting diodes^), and microelectronics (high electron mobility transis- 
tors^). Their utility in optoelectronics stems in large part from their wide range of direct 
band gaps - 0.7 eV in InN, 3.4 eV in GaN, and 6.3 eV in A1N - combined with their ability to 
form complete solid solutions, which provides the entire spectrum of emission wavelengths 
from 1800 nm (infra-red) to 200 nm (ultra-violet). In addition, the macroscopic polarization 
associated with the wurtzite crystal structure plays an important role in determining the 
electrical and optical properties. This is particularly important in heterostructures where a 
difference in polarization between layers induces an electric field at the interface, which in 
turn is screened by the formation of a two-dimensional electron gas (2DEG). The electron 
mobility in such polarization-induced 2DEGs can be much higher than that in traditional 
impurity-doped systems, in which the electrons scatter from ionized impurities. Thus high 
speed transistors known as POLFETs (polarization field effect transistors) can be produced 3 . 

Presently, there is growing interest in ZnO-based materials as alternatives to the nitrides. 
Indeed many properties of ZnO are similar, if not superior to those of GaN. Perhaps most 
importantly, and in striking contrast to GaN, single crystals of ZnO can be grown readily^, 
and used as substrates for the growth of thin film devices. This facilitates the production of 
much higher quality films using homoepitaxy, and circumvents the problems associated with 
dislocation formation from epitaxial mismatch that plagues GaN growth. In addition, the 
exciton binding energy of ZnO (60 meV) is twice as large as that of GaN, leading to stronger 
electron-hole recombination and thus improved optical efficiency^. Finally, the band gap of 
ZnO (3.4 eV) is very close to that of GaN, and it can likewise be varied systematically 
by alloying with MgO&£ or CdO 8 . However, since MgO and CdO do not occur in the 
wurtzite structure, values of the electronic and elastic properties, crucial for the modeling 
of hetero junction device performance^, cannot be obtained experimentally. 

In this work we calculate the spontaneous polarizations, piezoelectric coefficients and elas- 
tic constants of wurtzite-structure MgO, ZnO and CdO. We calculate the optimized struc- 
tures of the three materials using the standard local density approximation (LDA) to density 
functional theory as implemented in the VASP packag e) 10 ' 11 . We then use our recently- 
developed pseudopotential self-interaction corrected (pseudo-SIC) implementation^ of the 
density functional formalism to calculate the electronic properties, because the LDA fails 
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to obtain a band gap for wurtzite-structure CdO, preventing the calculation of its polariza- 
tion. The polarizations are obtained using the widely-used Berry phase metho d 13 ' 14 , and the 
piezoelectric coefficients and elastic constants are obtained from the calculated dependence 
of polarization and energy on the appropriate strain. Our polarizations are reported relative 
to the conventional reference structure of the zincblende structure with the same ratio of c 
lattice constant to unit cell volume^, and care is taken to avoid problems associated with 
the branch dependence of the polarization by extracting the so-called proper piezoelectric 
response for the e%i coefficients^. Our main result is that the polarization magnitudes and 
polarization gradients are all larger in the ZnO-based system than in the GaN-based sys- 
tem, suggesting ZnO as a promising candidate for the development of polarization field-effect 
transistors. 

Before reporting our polarization results, we compare our pseudo-SIC band structures 
with those that we obtain using the LDA (Figure [TJ, at the calculated LDA lattice pa- 
rameters. It is clear that our calculated pseudo-SIC band gap of ZnO (3.8 eV) is in good 
agreement with the experimental value (3.4 eVli), whereas the LDA shows the usual under- 
estimation (0.78 eV). CdO and MgO occur naturally in the rock-salt structure with band 
gaps of ~ 1 eV (indirect) or 2.4 eV (direct) for CdO 17 and ~ 7.7 eV for MgO 17 ; although 
our calculated band gaps for the wurtzite structure should not match exactly, they are of a 
similar magnitude. 




FIG. 1: Band structures of MgO, ZnO and CdO within (a) pseudo-SIC and (b) LDA approaches, 
at the calculated optimized LDA lattice parameters. Notice that the LDA severely underestimates 
the band gap, and results in a metallic ground state for CdO; this prevents the calculation of the 
polarization and piezoelectric coefficients. 
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P e 33 


e 3 i 




MgO 


-0.173 2.26 


-0.38 


II-VIs 


ZnO 


-0.048 1.34 


-0.57 




CdO 


-0.099 1.67 


-0.48 




A1N 


-0.081 1.46 


-0.60 


III-Vs 


GaN 


-0.029 0.73 


-0.49 




InN 


-0.032 0.97 


-0.57 



TABLE I: Polarizations, P, and piezoelectric constants, e 33 and e 3 i, in C/m 2 , for II- VI oxides 



(this work), and III-V nitrides (from Ref. 
lattice parameters. 



151 ). calculated at the theoretical LDA equilibrium 



When a wurtzite-structure semiconductor is constrained by epitaxial matching to a sub- 
strate or in a heterostructure, its total polarization is the sum of its intrinsic spontaneous 
polarization, P sp (that is the polarization that it would have in an unstrained bulk sample) 
plus the polarization induced as a result of the strain, P pz (e). The strain-induced component 
depends strongly on the strain, e, and hence on the lattice mismatch between the epitaxial 
layers; in the linear regime it is related to the piezoelectric tensor, e, by 
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Both the intrinsic and strain-induced contributions can be calculated accurately from first 
principles calculation a 13 ' 14 ' 15 ' 18 , and have been shown to be large in the III-V nitrides^. 

We begin by calculating the polarization and piezoelectric constants (e 33 and e 3i ) at the 
calculated LDA equilibrium structures for MgO (a = 3.26 A, c/a = 1.51, u = 0.398), ZnO 
(a = 3.20 A, c/a = 1.61, u = 0.378) and CdO (a = 3.60, c/a = 1.55, u = 0.391); our results 
are reported in Table HI For comparison we also report the literature values calculated for the 
nitrides at their theoretical LDA lattice parameters and u values^. (Note that our reported 
polarization for LDA ZnO (-0.048 C/m 2 ) differs slightly from the literature value (-0.057 
C/m 2 )^ because we use our pseudo-SIC approach; if we repeat our calculations using the 
LDA we reproduce the earlier value. Likewise, our pseudo-SIC piezoelectric constants are 
slightly different from literature LDA value o 15 i 20 .) 

Notice that the polarizations are larger in every case for the oxide materials than for their 
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nitride analogues. Perhaps more importantly, the polarization gradients between the end- 
point materials are considerably larger in the ZnO-based system. This is illustrated in Fig. |21 
where we compare our polarizations for MgO, ZnO and CdO with the literature values for 
A1N, GaN and InN, all at the theoretical equilibrium structures. The straight lines indicate 
the Vegard's Law projections for the polarizations of intermediate alloy compositions. 




FIG. 2: Calculated absolute polarization values of the II- VI oxides ZnO, MgO and CdO (this work) 



and their III-V nitride analogues GaN, A1N and InN (from Ref. 



In principle, the polarization in any other strain state can then be obtained from the 
equilibrium values shown in Table |l] using Eqn. ^ However, since the experimental bond 
lengths depend sensitively on the details of the structure, and could in fact be distorted 
beyond the linear regime described by Eqn. we also calculate the polarization from first 
principles in two other limits: First, at the experimental ZnO lattice parameters and atomic 
positions (a = 3.25 A, c/a — 1.60, u = 0.383^). This represents the situation in which only 
small amounts of Mg or Cd are doped into ZnO so that the structure does not change signif- 
icantly. And second, we constrain the in-plane lattice parameter a to that of experimental 
ZnO, and allow the structure to relax perpendicular to the plane (along the c direction) to 
simulate the situation in epitaxial films. 

Our results are shown in Table |HJ It is clear from the first column (calculated using the 
ZnO experimental coordinates) that, when the ions are constrained to have the same ionic 
positions, the polarizations of the three materials are rather similar. Therefore differences 
in polarizations in real systems must result primarily from the relaxations of the ions. The 
second column of data in Table ITT1 indicates that when MgO is allowed to relax out of plane, 
there is a large increase (almost double) in the polarization; this is caused by a contraction 
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p 




(i) (n) 


MgO 


-0.067 -0.152 


ZnO 


-0.073 -0.081 


CdO 


-0.073 -0.013 



TABLE II: Polarization values for the II- VI oxides calculated (i) at the ZnO experimental atomic 
positions and lattice parameters, and (ii) with the a parameter constrained to the experimental 
ZnO value, and the theoretical relaxed c parameters and atomic positions. 

of the c parameter during the structural optimization (c/a reduces from 1.60 to 1.52 and u 
changes from 0.382 to 0.395). Similarly the the large (factor of six) decrease which we find 
for CdO is caused by an out-of-plane expansion (c/a= 1.87, u = 0.357) to accommodate 
the applied in-plane compression. However the polarization of ZnO is almost unchanged 
from that at the experimental coordinates, because there is negligible change in the c lattice 
parameter on relaxation. Note also that the polarization of MgO is similar to that at its 
calculated equilibrium parameters (Table HJ), because the calculated equilibrium parameters 
are not dissimilar to those with the in-plane constraint. 

Finally, (Table ILTT]) we report the elastic constants (Cn, Ci 2 , C 33 , Ci 3 , and C 44 ) obtained 
by calculating the change in energy with strain for five different strain configurations. The 
details of the strain c onfig urations used and the procedure for extracting the elastic constants 
are described in Ref. |2l|. For comparison, the experimental values for ZnO, obtained using 
the Brillouin light scattering technique on polycrystalline films^ are also included in the 
table. Note that the agreement between our calculated values and the measured values for 
ZnO is very good. We expect our calculated values for hypothetical wurtzite-structure MgO 
and CdO to be similarly appropriate. 

In summary, we have calculated the spontaneous polarizations, piezoelectric coefficients 
and elastic constants for wurtzite-structure MgO, ZnO and CdO. For ZnO our results agree 
well with previous calculations and with experimental data; for MgO and CdO (which are 
experimentally inaccessible in the wurtzite structure) we provide predictions. We find that 
the absolute polarizations and the polarization gradients between end-point compounds are 
larger in the ZnO system than in the GaN system. 
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Cn C12 C33 


Cl3 


C44 




MgO 222 90 109 


58 


105 


calculated 


ZnO 217 117 225 


121 


50 




CdO 150 108 105 


61 


47 


measured 


ZnO 206 118 211 


118 


44 



TABLE III: Elastic constants for wurtzite-structure MgO, ZnO and CdO calculated at the theoret- 
ical equilibrium lattice parameters and atomic positions (this work) and measured using Brillouin 
light scattering (Ref. 
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